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ABSTRACT 

We investigate the effect of non-evaporating primordial black holes (PBHs) on the ionization and 
thermal history of the universe. X-rays emitted by gas accretion onto PBHs modify the cosmic re- 
combination history, producing measurable effects on the spectrum and anisotropies of the Cosmic 
Microwave Background (CMB). Using the third-year WMAP data and FIRAS data we improve ex- 
isting upper limits on the abundance of PBHs with masses > 0.1 Mq by several orders of magnitude. 

Fitting WMAP3 data with cosmological models that do not allow for non-standard recombination 
histories, as produced by PBHs or other early energy sources, may lead to an underestimate of the 
best-fit values of the amplitude of linear density fluctuations (as) and the scalar spectral index {us). 
Cosmological parameter estimates are affected because models with PBHs allow for larger values of 
the Thomson scattering optical depth, whose correlation with other parameters may not be correctly 
taken into account when PBHs are ignored. Values of ~ 0.2, ris ~ 1 and cts ~ 0.9 are allowed at 
95% CF. This result that may relieve recent tension between WMAP3 data and clusters data on the 
value of (78 ■ 

PBHs may increase the primordial molecular hydrogen abundance by up to two orders of magnitude, 
this promoting cooling and star formation. The suppression of galaxy formation due to X-ray heating 
is negligible for models consistent with the CMB data. Thus, the formation rate of the first galaxies 
and stars would be enhanced by a population of PBHs. 

Subject headings: cosmology: theory — cosmology: observations — early universe — cosmic microwave 
background — cosmological parameters — black hole physics 



1. INTRODUCTION 

During the radiation era, before the formation of 
the first stars and galaxies, small perturbations of 
the energy-density of the universe on scales compara- 
ble to the particle horizon may become gravitation- 
ally unstable. The outcome of the collapse can be 
the direct formation of primordia l black holes (PBHs) 
jHawking '1971'; 'Carr & Hawkingl Il974al: iMusco et all 
|2005; Harada & Carr 2005). If these PBHs form at a 
sufficiently high mass, they do not evaporate but instead 
begin to grow by accretion, producing x-rays. 

A population of sufficiently massive and numerous 
PBHs may provide an important and observable source of 
energy injection into the cosmic gas before the formation 
of non-linear large scale structures and galaxies. Sev- 
eral previous papers have addressed this scenari o 
[T98ltlGnedin et al.lfl995tlMiller fc Ostrikei]l200lD . Here, 
we improve on previous calculations by modelling the 
accumulation of dark matter around PBHs, the proper 
motion of PBHs and feedback effects including Compton 
drag. We simulate the ionization and thermal history of 
pr imordial plasma us i ng se mi-analytical calculations as 
in lRicotti fc Ostrikeil ()2004l ). The early energy injection 
by PBHs n iay produce observabl e distortions of the CMB 
spectrum ( |Battistelli et al.l 12000) and may affect CMB 
anisotropies. We fit cosmological models which include 
the effect of PBHs on the cosmic recombination history 
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to the 3rd year WMAP data (jSpergel et al.ll2006l ). This 
is done by modify ing the Monte-Carlo code COSMOMC 
(|Lewis fc Bridlell2002i) appropriately. Our main goal is 
to improve existing upper limits on the mass and abun- 
dance of non-evaporating PBHs. In addition, we find 
that the existence of PBHs may affect the best fit es- 
timates of cosmological parameters. In particular, the 
value of Te, (78 and Us can be underestimated in models 
which do not account for the existence of PBHs. This 
result is more general than the case of PBHs discussed 
in this paper. Any source of energy injection at early 
time that modifies the recombination history may lead 
to underestimating Te, Us and if the effect is excluded 
a priori when fitting the WMAP3 data. This scenario 
may ease recent tensions between the WMAP3 analysis 
that favors low values of (Tg ~ 0.74 (assuming a standard 
recombination history) and clusters data that may favor 
large r values of as ^ 0.9 (|Evrard et al.l 120071 ) . but see 
also (|Bode et al.ll2007t) . 

Finally, ionization from the X-rays of accreting PBHs 
will increase the amount of H2 significantly and the re- 
sultant extra cooling will enhance early star formation. 

This paper is organized as follows. In §[2] we review the 
properties of PBHs and the current constraints. In § [3] 
we introduce the basic equations for gas and dark matter 
accretion and for the accretion luminosity of PBHs. In 
§ |4] we evaluate the effect of local and global feedback 
processes. In § [5] we present the results of calculations 
of the ionization and thermal history of the IGM and 
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in § [6l data from WMAP3 and COBE are used to con- 
strain the mass and abundance of PBHs. A summary 
and discussion is presented in § [T] Throughout the rest 
of the paper we use the following cosmological parame- 
ters {h = 0.73, = 0.238, »b = 0. 0418,(78 0.74, = 
0.95) from WMAP3 (.Soergel et al.l[200fil 

2. PRIMORDIAL BLACK HOLES 

In the Newtonian regime, the theory of PBH formation 
can be understood in simple terms. The Jeans length 
in a static and homogeneous fluid with sound speed Vg 
is Rj ^ {c/vg)Rsch, where Rsch is the Schwarzschild 
radius of a black hole with density equal to the mean 
cosmic value. During the radiation era w c/\/3, thus 
Rj — > Rsch- This means that the self gravitating regime 
appears when the perturbation is very close to the black 
hole regime. The critical overdensity needed to trig- 
ger PBH collapse is w ;S 1/3, where P = wpc^ is th e 
cosmic equation of state ()Carijll975t iGreen et al.ll2004D . 
The typical mass of a PBH is approximately equal to or 
smaller than the mass within the particle horizon at the 
redshift of its formation, zf. 

Mpth J. f 1 + Zeg\'^ ^ f P '^'^ 



(1) 



MniZeq) \1 + Zf J \ Qpbh , 

Here, Jhot — -^^pfch/-M^ff (z/) < 1 is the ratio of the 
mass of the PBH to the mass of the particle horizon at 
z = z/, /3 is the density parameter of PBHs at z = Zf, 
Mnizeq) ~ 3.1 X 10^^ Mq is the mass of the horizon 
at the redshift of m atter-radiation equa lity Zeq ~ 3000 
(for ^ 0.127. ISpergel et al.|[200i) . The relation- 

ship /3(1 -|- Zf) = Qpbhi^ + Zeq), valid for non-evaporating 
PBHs, relates the density parameter of PBHs today, 
^pbh{M) ~ Ppbh/Pcrit, to the ouc at the time of for- 
mation, (3{M). For example, if a tiny fraction (3 ~ 10~^ 
of the cosmic energy- density collapses into PBHs during 
the quark-hadron phase transition at t 10~^ s, it fol- 
lows from equation ^ that about 30% of the present day 
dark matter is made of PBHs with mass Mpth ^ 1 Mq 
(we assumed fnor = !)• Nothing prevents the dark mat- 
ter from being a mixture of weakly interacting particles 
(WIMPs) and PBHs with fldm — ^pbh + ^wimp- How- 
ever, equation ([1]) shows that fine tuning of the value of 
P is required in order to have flpbh ~ ^wimp- Thus, it 
is often considered more likely that dark matter is dom- 
inated by either WIMPs or PBHs. PBHs may have an 
extended mass range if they form at different epochs and, 
even if they all form at one epoch, simul ations show that 
their mass distribution is broad {e.g., IChoptuifl 119931 : 
lEvans &: Colemanlll994t ). Thus, present-day dark mat- 
ter may be composed of a mixture of relatively massive 
and tiny PBHs, some of which may be completely evap- 
orated or have left Planck-mass relics. In practice, using 
astrophysical tests only, sufficiently small PBHs would 
be virtually indistinguishable from WIMPs. It is worth 
pointing out that no fine tuning of (3{zf) is required to 
produce ^Idm = ^pbh 1- Any value of /3 < 1 and 
Zf is allowed but, of course, different values of the en- 
tropy of the Universe and different redshift of equiva- 
lence would be produced. Specifically, the radiation den- 
sity parameter today is related to /3 and z/ by the re- 
lationship Qradfi = ^dm/{l + Zeq) ^ /3""^(1 + Zf)"''' . In 

addition, a scenario in which PBHs produced after in- 
flation evaporate leaving only Planck mass relics may 



explain why flradflh^ ~ 4.35 x 10 ^ and provide a justi- 
fication for the large value of the entropy of the Universe 
() Alexander fc Meszaroil2007t) . 

PBHs are a unique probe of the early universe, of high 
energy processes and of quantum gravity. The small mass 
scales at which PBHs may form are inaccessible to Cos- 
mic Microwave Background (CMB) experiments. As of 
today there is no solid evidence for the existence of PBHs, 
but their presence would be very difficult to detect even if 
they constitute the bulk of the da rk matter. Ea rly results 
from the MACHO collaboration (i Alcock et allfeOOO l sug- 
gested a possible detection of solar mass objects consti- 
tuting about 20% of the Galactic halo mass. Recent re- 
sults seem to disfavor this claim {e.g., iHamadache et ahl 
l2006f) . The goal of the present work is to better con- 
strain the abundance of relatively massive {M 0.1 Mq) 
PBHs by modeling their effect on the ionization history 
of the universe. To further motivate the present study 
we briefly review the current status of theoretical works 
and observational limits on the e xistence of PBHs. For 
a more comprehensive review see [C!anl(l2005h . 

2.1. Current observational limits on PBHs 

Many physical processes may lead to the formation 
of PBHs. For example, PBHs may form from pertur- 
bations after inflation, during phase-transitions of the 
cosmic equation of stat e or from t opological defects 
(|Carr fc Hawkind Il974bl: ICaTJ [20051 e.g.). Theoreti- 
cally, it is unclear what is the largest allowed mass of 
PBHs. Inflationary theories predict an almost scale- 
invariant {n ~ 1) initial spectrum of perturbations. 
This spectrum has the remarkable property that all 
perturbations have the same amplitude when they en- 
ter the horizon. Thus, in this case, it may appear 
that the probability of PBH formation is almost inde- 
pendent of their mass. Small-mass PBHs would have 
a larger probability of formation if the spectrum were 
slightly blue (with n > 1) and vice versa if the spec- 
trum were red, as see ms to be indicated by availab le data 
JSpergel et al.' '20061 [Chongchit nan fc EfstathToul 120071) . 
IJedamzik (199:^ have shown that during a first-order 
phase transition the cosmic equation of state may be- 
come softer {w <C 1/3). As a result, the formation 
of PBHs may be substantially enhanced at that par- 
ticular mass scale. PBHs with masses ~ 1 Mq may 
form during the QCD (quark-hadron) phase transition 
at t ^ 10^^ s, or PBHs with mass 10^ Mq may form 
during — annihilation era. PBHs with masses of 
100 — 1000 Mq may be produced in two-stage inflation- 
ary models designed to fit the low WMAP quadrupole 
(jKawasaki ct al. 2006). One interesting feature of this 
model is a huge bump in the power spectrum at kpc 
scales with overdensity 6 > 1. Other possible scenar- 
ios for PBH formation involve non-Gaussian perturba- 
tions produced by cosmic strings, cosmic string collapse 
or bubble co llisions following second-orde r phase tran- 
sitions (e.g.. iPolnarey fc Zem bowicz 1991|; iRubin et al.l 
I2001t IStojkovic et al.l l2005l : iNozari .2007i) . It is gener- 
ally considered unlikely that PBH formation occurred 
after t = 1 s, when Mpbh ~ Mh ^ 10^ Mq, because 
the physics in this domain is sufficiently understood and 
their formation would affect primordial nucleosynthesis. 

Observationally, the abundance of PBHs is well con- 
strained only at very small masses. PBHs with masses 
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smaller than 5 x 10 g should not exist today because 
they evaporate in le ss than a Hubb le time by emission of 
Hawking radiation (lHawkindll974f ). The radiation from 
evaporating PBHs affects nucleosynthesis and the CMB 
spectrum and may overproduce the observed gamma-ray 
background. Hence, upper limits on the abundance of 
PBHs with masses 1 g < Mp^/j < 10^^ g are quite strin- 
gent: f3{M) - 10-20_ 10-22 fe.q. [Caril[2003l ). The num- 
ber density of PBHs with masses larger than 10^^ g is 
poorly constrained because they emit negligible amounts 
of Hawking radiation (their timescale for evaporation is 
longer than the age of the Un i verse today). Dynamical 
constraints ("Lacev & Ostriked 119851 : iCarij 119941 : iMoord 
[l993; Jin ct al. 2005), effects on the matter power spec- 
trum (lAfshordi et al.l 120031) and statistical s tudies of bi- 
nary stars in the Galactic halo (|Yoo et al.|[2b04) . have 
been used to rule out PBHs more massive th an 1000 Mp, 
as ma in constituent of the dark matter, but (|Mack et al.l 
120071 hereafter Paper I), summarizing these constraints, 
find that a domain remains for which massive PBHs 
can make a significant contribution to the dark mat- 
ter. Searches for microlensing events toward the Large 
Magellanic Cloud by the MACHO and EROS collabo- 
rations so far provide the most stringent constraints on 
the existence of non-evaporating PBHs. They have ruled 
out PBHs as the bulk of the galactic dark matter in the 
mass r ange 10"'^ Mp < M < 30 M© ()Alcock et al.lll998l 
1200 If) . lAlcock et al.l ()2000l ) have claimed a positive detec- 
tion of Massive Compact Halo Objects (MACHOs) with 
M ~ 0.1 - 1 Mq constituting about 20% of the Milky- 
Way dark matter halo. More recent works have not con- 
firmed this result (jHamadache et ani2006[ ) but this mass 
range is of particular interest because it is about the mass 
of PBHs that may for m during quark-h adron phase tran- 
sition as proposed by f Jed amziMll997f ). The gray curves 
in Figure [9ljleft) summarize the current status of obser- 
vational upper limits on the abundance of PBHs with 
Mpth > 10^^ g. The two black curves labeled "FIRAS" 
and "WMAP3" show the new upper limits derived in the 
present work. 

2.2. PBH growth and "clothing" dark halo 

PBHs more massive than 10^ Mq may only exist 
if substantial accretion takes place after their forma- 
tion or if many PBHs which form in a cluster merge 
into a larger on e ([Mcszaros 1975a; Khlo pov et al.l [20051 : 
IChisholriil|2006[ ) . In the Newtonian approximation, PBHs 
can only grow if their mass is comparable to the hori- 
zon mass. In this case their mass increases at the 
same rate as the horizon mass (IZel'Dovich fc Novikovj 
Il967f). Relativistic calculations by iCarr fc Hawking 
( 1974bf ) have shown that growth by accretion is always 
slower than the growth of the horizon mass. Thus, it 
is generally accepted that the mass of PBHs does not 
increase significantly after their formation. In mod- 
els in which the e quation of state becomes stiff w ^ 
1 (|Lin et all [l97l or in brane world cosmologies in 
which extra dimensions are sufficiently large, growth 
by accretion could be significant and in principle may 
lead to the production of massive PBHs startin g from 
smalle r s eeds ([B ickncll fc Henrikscn 1978; Guedcn s et al.l 
I2002allbl : lMaimndai1l200a iTikhomirov fc Tsalkoull200l . 

Scenarios in which the bulk of the dark matter is not 
made of PBHs are particularly interesting because, after 



the redshift of matter-radiation equivalence, each PBH 
seeds the growth of a dark halo that, over time, becomes 
considerably more massive than the PBH in its center 
(Paper I). The halo gravitational potential may increase 
the gas accretion rate onto the central PBH by several 
orders of magnitude. If the accretion is near the Ed- 
dington rate, the radiation efficiency of PBHs may be 
large. At later times, when PBHs are accreted by the 
massive halos of galaxies, they likely lose much of their 
dark matter "clothing" and become "naked" again. 

3. BASIC EQUATIONS 

Here we review the theoretical framework for the ac- 
cretion of dark matter and baryons by primordial black 
holes. A point mass, M, immersed in an hydrogen gas 
with constant number density rigas and sound speed Cg 
that moves at velocity v with respect to the gas accretes 
at a rate 

A/b = X-iTTmHrigasVeffr^, 

where rs = GMv~J'^ is the Bondi-Hoyle radius and 

Veff = iy^ + clf/'^ l|Bondi fc Hovlelfl94l . The numeri- 
cal values of the Bondi radius and accretion rate are: 

M.»2xlO-g.-.A„,„(^)'(^JiiL^)l3) 

The mean cosmic gas density is 

n„as ^ 200 cmr'^ (^-^\ ■ (4) 

gas ^QQQ J V ) 

Before decoupling the gas temperature is roughly equal 
to the temperature of the CMB. Afterwards the gas tem- 
perature decreases adiabatically due to the Hubble ex- 
pansion of the Universe. The temperature of the IGM 
(neglecting for the moment the feedback produced by 
accreting PBHs) is well approximated by the fitting for- 
mula 

T„as = (2730 K) ( ^) ^ , (5) 

where a = l/(l-|-z) is the scale parameter, a^ec 
is the scale parameter at decoupling where z^ec — 
132(rJhft.V0-022)2/5 and (3 = 1.72. The gas sound speed 
is Cs = (5.7 km s~^){Tgas/'273Q)^^'^ . Thus, from equa- 
tion ([5]) we have 

/l + z\'/' 

«(5.7kms-i) f^^j forz » Zdec - 132. (6) 

For spherical accretion onto a point mass, assuming a 
non-viscous fluid, the parameter A is of order unity. In 
the present study, we need to take into account the effects 
of the growth of a dark halo around PBHs, the Hubble ex- 
pansion and the coupling of the CMB radiation to the gas 
through Compton scattering (Compton drag). All these 
effects can be folded into the calculation of the accretion 
eigenvalue A. In a companion paper, (Ricotti 20Q^ here- 
after Paper II), we have derived analytical relationships 
for A for the cases of a point mass and an extended dark 
halo, including the aforementioned cosmological effects. 
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Fig. 1. — (Left). Ratio of the baryonic to dark matter power spectrum as a funtion of wavenumber k. Each curve, from bottom to the top, 
corresponds to scale factors a = 0.001,0.002,0.005,0.01,0.05,0.1, respectively. (Right). The solid curve shows the mean relative velocity 
between the dark matter and baryons in a sphere of comoving radius rg . The dashed and dotted curves show the velocity dispersion within 
a sphere of comoving radius ro for the baryons and for the dark matter, respectively. 



In the following calculations we will use the results from 
Paper II to derive the accretion rate onto "naked" PBHs 
(appropriate only if PBHs constitute the bulk of the dark 
matter) and for the general case of accretion onto PBHs 
including the growth of a "clothing" dark halo of mass 

To obtain this result, we must first estimate the dark 
matter accretion (i 33.ip and the effects of PBH velocity 
ft j3.2p . the angular momentum ft ^3.3p of the dark and 
baryonic matter and PBH clustering f i|3.4[) . From the 
gas accretion rate, we can then estimate the accretion 
luminosity, which will ultimately be responsible for al- 
tering the evolution of the IGM. 

3.1. Growth of a dark halo around PBHs 

If PBHs do not constitute the bulk of the dark matter, 
they seed the accumulation of a dark halo of WIMPs 
which grows pro portionally to t^/^ (jBertschingeij [l985l : 
iMack et al.ll2007l ) . A similar result holds if the dark mat- 
ter is made only of PBHs but with a very broad distri- 
bution of masses. 

After the redshift of equivalence the mass of the dark 
halo surrounding PBHs grow proportionally to the cos- 
mic scale parameter: 



Mh{z) = (|)^Mpbh 



Z+l V 

1000 / 



(7) 



where 0i = 3 (see Paper I). The growth of the PBH 
during the radiation era is of order of unity. After z ~ 30, 
depending on the environment, the dark halo may stop 
growing: if a PBH evolves in isolation it can continue to 
grow, otherwise it will lose mass due to tidal forces as it 
is incorporated into a larger galactic dark halo. 

If Mh > Mpbfi and the accretion is not perfectly spher- 
ical, the dark halo has a self-s imilar power-law d ensity 
profile p cx r" with a ~ —2.25 (|Bertschingeijll985l) trun- 



cated at a halo radius 
ri, = 0.019 pc 



Mh 



1/3 



(8) 



1 MqJ \ 1000 

where is about one third of the turn-around radius 
(Paper I and II). 

3.2. Proper motions of PBHs 

The motion of the accreting mass relative to the sur- 
rounding material strongly affects the accretion rate (see 
equation [5]). This proper motion is determined by the 
relative amplitude of inhomogeneities of dark matter and 
gas. More precisely, assuming that PBHs behave like 
dark matter particles and neglecting binary interactions, 
we estimate their relative velocity with respect to the 
gas in the linear and non-linear regimes. We show that 
at redshifts z ^ 30 — 50, due to the imperfect coupling of 
gas and dark matter perturbations {i.e., Silk damping), 
the peculiar velocity of PBHs with respect of the gas is of 
the order of the gas sound speed. At redshifts z < 30 the 
growth of non-linear perturbations dominates the motion 
of PBHs. The velocity of PBHs falling into the potential 
wells of the first galaxies is sufficiently large to virtually 
stop the accretion of gas from the intergalactic medium 
onto PBHs, until they came to rest at the center of the 
halo into which they are accreted. 

3.2.1. Linear regime 

In linear theory, gas and dark matter perturbations 
are coupled. Before the redshift of "decoupling" at z ^ 
100, the growth of gas inhom ogeneities on small scales is 
suppressed by Silk damping (|Silklll968l ). At these scales 
and redshifts, the gas flow lags behind the dark matter 
until it is able to catch up at later times. This process 
is illustrated quantitatively in the left panel of Fig. [H 
which shows the ratio of the gas to dark matter power 
spectra in our fiducial ACDM model. The curves, from 
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bottom to the top, show the ratio of the power spectra 
of the baryons to dark matter calcul ated using the code 
"lingers" in the Graphic 1 package ijBertschingeil Il985l ) 
at scale parameters a = 0.001,0.005,0.01,0.05, and 0.1, 
respectively. 

^From the power spectra of dark matter and baryons 
we calculate the ensemble average of the center-of-mass 
velocity of a patch of t he universe of comoving radius ro 
(lOstriker fc Sutol[l990h : 
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27r2 



Pi{k)wl{k,a)wf{k,ro)dk, (9) 



where rim is the cosmological density parameter, Ws and 
wi are window functions (here we use "top hat" window 
functions) and a is a small scale smoothing of the pertur- 
bations. The choice of the value of a is not critical as long 
as a <C ^Q. The index i — dm, bm refers to dark matter 
and baryons, respectively. The ensemble average of the 
velocity variance within a patch of comoving radius is 
calculated in a similar fashion: 



P,{k)wl{k,a)[l - wf{k,ro)]dk. 



(10) 

This equation will be used later (§ 13. 3p to estimate the an- 
gular momentum of gas and dark matter accreting onto 
PBHs. The "cosmic Mach number," Mi = {Vi)/{(Ti), is 
independent of the uncerta inties on the nor malization of 
the power spectrum (Ostriker fc Sutol[l990f) . 

Since the flow of the gas and the dark matter trace 
each other, their mean relative velocity is (Vrei) = 
{Vdm) — {Vbar)- The right panel of Fig. [1] shows the 
mean relative velocity {Vrei) of PBHs with respect to 
the baryons (solid curves), and the velocity dispersion 
{o'dm,bm) for the gas (dashed curve) and the dark matter 
(dotted curve) within a sphere of comoving radius rg at 
z — 1000. Power-law fits for (a) as a function of ro and 
z are: 



Cfem ~ (0.35 km s 



(Tdm« (1-58 km s ^) 



, 0.85 / -, , \ -1 

ro \ f 1 + z 



1 Mpc J V 1000 

, 0.85 

''O 



1 Mpc 




fll) 



.(12) 



The fits are accurate to about 5% for rg < 1 Mpc in the 
redshift range 50 < z < 2000. 

In § 13.61 we discuss in more detail the radiative effi- 
ciency of accreting PBHs. We anticipate that the radia- 
tive efficiency e and the accretion luminosity L depend 
on the accretion rate m: L (x v'^^^ if to < 1 and L cx v~^j: 

if m 1, where the effective velocity Veff is defined in 
§ |3] (see equation [3]). The effect of the linear velocity 
field, which is Gaussian, on the mean accretion luminos- 
ity of PBH is given by the statistical ensemble average of 
v~"j: weighted by the distribution function of the PBHs 
relative velocities with respect to the gas: 



fM{v,(j)dv 

(c2 +„2)a/2^ 



(13) 



where a = 3 or a = 6, depending on the value of the 
accretion rate. The linear velocity field is Gaussian; thus 
the distribution of the moduli of the 3D velocity field is 
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Fig. 2. — Luminosity weighted effective velocity of PBHs (thick 
curves): (t)e//)A (sohd curve), is weighted assuming I oc rr? and 
(v^ff)g (dashed curve), is weighted assuming / oc m. See the text 
for details. The thin curves show the variance of the velocity dis- 
tribution, (Vrei) (dashed line) and the gas sound speed, Cs (dotted 
line), respectively. 

a Maxwellian, /Af(w,cr) with a = (Vrei)- Finally, inte- 
grating equation (fT3|) and defining Mpbh — {Vrei)/cs, we 
find 



{Veff)A 



for a = 6 and 




for Mpbh > 1, 
for Mpbh < 1, 




f In(fA^pfc^) 



for Mpbh > 1, 
for Mpbh < 1, 
. (15) 

for a = 3. The calculations presented in this section on 
the effect of linear perturbations on the accretion lumi- 
nosity of PBHs are summarized in Fig. [51 The thick solid 
and dashed curves show {vef f)B from equation (fTS)) , and 
{veff)A from equation (fT4|) respectively, as a function 
of redshift. For comparison, the thin dashed and dot- 
ted curves show (Vrei) and the sound speed of the gas, 
respectively. Neglecting feedback effects, the PBH ve- 
locity is comparable to the gas sound speed at redshifts 
z ^ 500. On average, the motion of PBHs in the linear 
regime reduces the gas accretion rate by a factor of a few 
with respect to the static case. 

3.2.2. Non-linear regime 

The velocity of PBHs falling into the gravitational 
potential of nearby galactic halos can be roughly esti- 
mated assuming that it is comparable to the circular 
velocity of virialized halos of 2a density perturbations: 
Vp ~ Vc{2cr,z). We c alculate Vc(2a, z) using the Press- 
Schechter formalism (jPress fc Schechteri[l97l . Adopt- 
ing a top-hat window function to calculate the variance of 
density perturbations, we find cr(M) « 10.2-0.79 log(M) 
for masses M < 10^° M0. Hence, for the mass range of 
interest. 



8 X 10^^ M0)exp[-1.8(z-H 1)], 



(16) 
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is a sufficiently accurate approximation for the mass of 
2(7 perturbations as a function of redshift. From equa- 
tion (jl6p we calculate the circular velocity and thus typi- 
cal proper velocity of PBHs induced by non-linear struc- 
tures, as a function of redshift: 



(17 km s 



-'if- 



2a 



Z + 1 



(17) 



The fraction of the dark matter and hence of PBHs 
in virialized halos with mass > M„i„ is fviriz) = 
1 — erf (i^min 7^/2) , where erf is the error function and 
u = 5c/ cr{Mmin, z). For a fraction of PBHs we as- 
sume that they have a velocity Vp relative to the gas and 
that the gas has overdensity 5 — 200. The remaining 
PBHs are in the intergalactic medium and are subject to 
global and local thermal feedback (see § 2]) . 

Dynamical friction may allow the most massive PBHs 
to spiral to the center of the halos in less than a Hub- 
ble time (see § [7]) , where they meet favorable conditions 
for accretion. For smaller PBHs inside virialized halos, 
their contribution to the accretion luminosity is impor- 
tant only in the redshift range \Q ^ z ^ 30. At z > 30 
the fraction of PBHs in virialized halos is small. At 
z < 10, the first massive halos form and the rate of 
gas accretion onto PBHs decreases rapidly due to the 
increasing relative velocity of the PBHs with respect to 
the gas: Vc ;^ 10 kms~^. During this redshift interval 
only a fraction ~ 10 — 20% of PBHs is inside virialized 
halos. Most PBHs at z ~ 10 are still in the intergalactic 
medium. The gas accretion onto the PBHs in the inter- 
galactic medium is strongly suppressed at z < 10 due to 
global thermal feedback (see Fig. [B]). 

3.3. Angular momentum of accreted material 

The angular momentum of the accreting gas deter- 
mines whether a disk forms or the accretion is quasi- 
spherical. In turn, the geometry of gas accretion deter- 
mines the radiative efficiency e(m), i.e., the efficiency 
of conversion of gravitational energy into radiation (see 
§ 13. 6p . The angular momentum of accreting dark mat- 
ter determines the density profile of the dark halo en- 
veloping a PBH and whether the mass of a PBH can 
grow substantially by accreting a fraction of its en- 
veloping dark halo. The angular momentum of the 
gas and dark matter accreting onto PBHs can be esti- 
mated from equation (fT^ for the mean values of the 
velocity dispersion within a comoving volume of radius 
ro (see right panel of Fig. [1]). For the gas we esti- 
mate Ubmir B ,com) sX the Bondi radius (see equation [3]), 
and for the dark matter we estimate <Jd7n{rh,com) at the 
the turn-around radius (see equation [5]). If we ne- 
glect the proper motion of PBHs, the Bondi radius of 
a "naked" PBH is approximately constant at redshifts 
z < 100: rs,,o,„ - (1.3 x 10"^ Mpc)(Mp6,,/lMQ). The 
proper motion of PBHs produces a small reduction of 
the Bondi radius that we parameterize with the function 
^(z) = max[Mpbhiz), 1], where Mpth = {'"eff)/cs, is the 
PBH Mach number. Replacing Mh with the mass of the 
PBH in the previous expression, we obtain an estimate 
of the Bondi radius including the enveloping dark halo. 
Finally, from equation (|12p we obtain 

■ ■ ' ■ Mh ^ 



O'bv 



' 0'f)m,oC(^) 



-1.7 



1 + z 
1000 



1 Mr. 



, (18) 
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Mh 



0.28 



(19) 



1000 ) V 1 , 

where Ohm.Q = 3.8 x 10^^ km s^^ and adm.o — 1-4 x 
10-* km s \ 

Let's first consider the accretion geometry of the gas. 
Applying conservation of angular momentum we find 
that the rotational {i.e., tangential) velocity of the gas at 
a distance r from the black hole is v{r)r ~ crhmrB, where 
rB is the Bondi radius. Neglecting relativistic effects, if 
the velocity is much smaller than the Keplerian velocity 
in the proximity of the black hole, then the accretion is 
quasi-spherical; if vice versa, a disk can form. The Ke- 
plerian velocity at radius r is Vkep = c{Rsch/'fy^'^ < c, 
where Rsch = 2GM/c^ is the Schwarschild radius of the 
PBH {Rsch ~ 3 km for a 1 Mq black hole). Thus, the gas 
accretion is quasi-spherical if abm < '2D£^{z)'^c'j./c. We 
have assumed Cs oc {1 + zy^^ and the constant 13 ^ 1 — 10 
takes into account relativistic corrections. We conclude 
that the angular momentum of the gas is negligible when 
the halo mass is 



Mh < il7A6MQ)D'-'^^{z) 



^4.33 



1 + 2 

1000 



2.35 



(20) 



Using A4pbh{z) calculated in § 13.2.11 and D ~ 10, the 
flow is quasi-spherical for any halo mass Mh ^ (1700 — 
17000) M0, almost independently of redshift. At red- 
shifts z < 100, the dark halo mass is about 30 times 
the mass of the PBH in its center. We conclude that an 
accretion disk forms at redshifts z < 100 around PBHs 
more massive than Mpbh ~ 500 Mq. Thermal feedback 
effects, which are important at z < 100, will increase our 
simple estimate of the critical mass for disk formation. 
So far we did not take into account feedback effects and 
Compton drag which are important in reducing the gas 
angular momentum. Previous studies have found that 
due to Compton drag the flow is nea rly spherical at r ed- 
shifts z > 100 {e.g., Loeb 1993; Ume mura et al.lll993l ). 

Similar arguments apply to the dark matter compo- 
nent. The flow is sufficiently spherical for the dark 
matter to be directly accreted onto the PBH if adm < 
cDRsch/fh- This inequality is satisfied for PBHs with 



Mpbh > (7.18 X 10^ Mq)D-^-' (^i±^^ 



-2.9 



M, 



pbh 



13\ n-2.6 



M 



H 



>(7.76 X 10"*'^)i:> 



l + z \ 
1000 / 



(21) 



in terms of the Horizon mass Mh{z). Including relativis- 
tic corrections {D ~ 10), only PBHs with mass ^ 10'^ 
Mq a,t z ^ 1000 (note the steep dependence of the crit- 
ical mass on redshift) may grow substantially from di- 
rect accretion of dark matter. The geometry of the dark 
halo becomes more spherical with the increasing mass 
of the PBH, since the ratio of the halo radius to the 
Schwarzschild radius, rh/Rsw « M/^^^j'^ , decreases with 
increasing mass. We conclude that direct accretion of 
dark matter into PBHs is typically negligible. In most 
cases the density profile of the dark halo is well described 
by the Bertschinger self-similar solution with log slope 
a = 2.25. 
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Fig. 3. — The dimensionless accretion rate of baryonic matter 
redshift. The curves from bottom to top refers to Mpi,i^ = 1, 10, 100 
Xe = 10~^; the right panel forxe = 1. Here, thermal feedback and 
motion of the PBH due to linear density perturbations is included. 

3.4. PBHs Clustering 

Recently it has been s uggested by several authors {e.g., 
iDokuchaev eTan[2?)0l [C^[200l [C^isholnj[200l that 
PBHs have high probability of forming clusters or bi- 
naries. A precise calculation of the effect of clustering is 
complex and beyond the scope of the present paper. Nev- 
ertheless, there are simple arguments that we can use to 
show that clustering may increase the radiative efficiency 
of PBHs. 

We consider two regimes determined by the relation- 
ship between the typical size of the cluster (or the the bi- 
nary separation), rd, and the Bondi radius of the whole 
system, rB,tot- 

1) If rd > rB,tot, the orbital velocities of the PBHs are 
smaller than the effective translational velocity of the 
system. Hence, we can neglect the orbital velocities to 
estimate the accretion rate. To a first approximation in 
this case, we can ignore the fact that the PBHs are clus- 
tered or in binary systems for the purpose of calculating 
their accretion luminosity. 

2) If rd < rB,tot, the orbital velocities of the PBHs are 
larger than the effective translational velocity of the sys- 
tem. Thus, the material is not accreted directly onto the 
PBHs but near the center of mass of the system and sub- 
sequently onto the PBHs. Because the PBHs are orbiting 
the center of mass of the cluster, the angular momentum 
of the accreted material is large and so the accretion ge- 
ometry is disk-like rather then spherical. The formation 
of a disk can increase the radiative efficiency by a fac- 
tor of ten (see Further work is needed to better 
understand the accretion rate of PBHs in this regime. 

3.5. Gas accretion rate 

We now have all that is needed to estimate the ac- 
cretion rate onto PBHs from equation ([3]). We define 
the dimensionless accretion luminosity / = Ltoi/LEd, 
where LecI = 1-3 x 10^^{Mpbh/l Mq) erg s~^. The ra- 



onto a "naked" PBH (without enveloping dark halo) as a function of 
300, 10^, 10*, lO'^ M0. The left panel is for a gas with electron fraction 
the contribution from PBHs inside virialized halos are neglected. The 

diative efficiency e determines the accretion luminosity 
for a given accretion rate: L^oi — eMbC^. Thus, defin- 
ing the Eddington accretion rate as MecI = Leci/c^ = 
1.44 X 10^"^ {Mpbh/MQ) g s~^ and the dimensionless ac- 
cretion rate as m = Mb/MEd we have I — em. 

The dimensionless Bondi-Hoyle accretion rate for a 
PBH without a dark halo is: 



m = (1.8xlO-^A) 



l + z V 
1000 / 



M, 



pbh 



IMr, 



Veff 



5.74 km s 



(22) 

This equation is valid only if PBH constitute all the 
dark matter: fp^h = ^pbh/^dm — ^- If the accret- 
ing object is a point mass (i.e., a "naked" PBH), A de- 
pends only on the value of the dimensionless gas viscosity 
j3 = p^ffrB/cs-. if /3 < 1 we find A - 1 and if /? > 1 we 

have A ~ 1//3. The fit to the accretion eigenvalue for an 
isothermal gas is given by the formula (see Paper II), 



exp 



9/2 



3 + 13 



0.75 



(23) 



where the dimensionless sonic radius {xcr = rrc/fB) is 
Xcr = [—1 + (1 + /3)^^^]//3- The effective gas viscosity, 
P'^ff = (3 + H, is the sum of two terms. The first term, 
due to Compton drag, is /3(z) = 2.06 x 10~^^a;e(l -|- z)^ 
s~^, and the second, due to the cosmic expansion, is the 
Hubble parameter H{z). It is easier to understand the 
physical meaning of the dimensionless viscosity due to 
the Hubble expansion in terms of the recession velocity 
of the gas due to the Hubble fiow at the Bondi radius: 
vh ~ rBH{z). As expected, if vh/cs = ^Hubble < 1 
we can neglect the cosmic expansion. The dimensionless 
viscosity depends on the redshift, the mass of the PBH 
and the temperature and ionization fraction of the cosmic 
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Fig. 4. — Same as in Fig. |3]but including the growth of the dark halo surrounding the PBH {a 
bottom to top refers to Mpi,^ = 1, 10, 100, 300, 10^, 10*, 10^ Mr; 
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2.25 and 0i = 3). The curves from 
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5.74 km s-i 
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as a point mass in terms of accretion rate, sonic radius 
and dimensionless viscosity. The numerical value of x is 



X = 0.22 



(24) 



l + z \ 
1000 / 



1 kms 



(28) 



Fig. [3] shows the accretion rate m from equations 
((23|) as a function of redshift. The curves from bottom 
to top refer to Mpbh = 1, 10, 100, 10^, lO'^, 10^ Mq. We 
assume a gas with constant electron fraction Xe = 10^"^ 
after recombination in the left panel and Xe = 10^^ in 
the right panel. Thermal feedback is neglected and re- 
combination is instantaneous. 

If fpbh < 1 we need to consider the effect of a grow- 
ing dark matter halo with mass Mh{z) = <t>iMpi,h[{'^ + 
z)/1000]"^ For = 3 we find 



m = (0.016A) 



1000 / 



phh 



Veff 



5.74 km s 

(25) 

In the simulations presented in § [SJ we impose that the 
dark halo stops growing when all the available dark mat- 
ter has been accreted; i.e., when fp},h4'i[{^+z) — 
1. The gas accretion onto an extended dark halo of mass 
Mh and power-law density profile with logarithmic slope 
2 < a < 3 is related to the one for a point mass (with 
Mphh = Mh) by the following scaling relationships (see 
Paper II): 



(f ) r.r. (26) 



where x = tbI^k < 2, p = 3 — a, and 
T = f 1 + lO/jC')) ~' exp (2 - 



X)(f 



(27) 



When the Bondi radius tb is larger than twice the radius 
of the dark halo (x > 2), the dark halo behaves the same 



We adopt the parameters a = 2.25, 0; = 3 and given 
in equation ([5]), valid for a quasi-spherical dark halo in 
which the dark matter does not fall directly onto the 
PBH. Figure [4] shows the accretion rate, to, from equa- 
tions ((25 |l -([26 |) as a function of redshift, with the curves 
representing the same masses as Fig. [3] with the same 
simplifying assumptions going into the simulation. 

3.6. Accretion luminosity 

The radiative efficiency e of accreting black holes is 
known to depend on the accretion rate to and the geom- 
etry of the accretion flow. Observations and theoretical 
models suggest that the radiative efficiency approaches a 
constant value e '--^ 0.1 for accretion rates to 1. If the 
accretion rate is m 1 the radiative efficiency is smaller 
by a factor e cx m (e.g.. [Park fc Ostrikerj|2001f ). Depend- 
ing on the angular momentum of accreted material, a 
thick disk may form. Otherwise the accretion geometry 
is spherical. The later case is the most conservative as the 
accretion efficiency is minimum. In addition, this case is 
appropriate for most PBH masses which have to < 1 (see 
?! 13. 3p . We consider the following accretion regimes: 
Thin disk: If m > 1 we assume that the gas accreted by 
the PBH forms a thin accretion disk and has radiative 
efficiency e = 0.1. We assume that the bolometric lumi- 
nosity cannot exceed the Eddington limit and that, due 
to complex feedback effects that we do not attempt to 
model here, the PBHs accrete near the Eddington limit 
only for a fraction of time fduty Observations of AGN at 
redshifts z < 3 show that ty pically fduty is about three 
percent (jShaplev et al.l 12003. e.g.). However, for PBHs, 
the duty cycle can be substantially different as the nature 
of the feedback mechanisms which regulate the cycle are 
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not well understood. We let fduty be a free parameter 
and we show that the limits on the PBH abundance scale 
linearly with it. In summary, 

I = fpbh X min(0.1m, 1) if m > 1. (29) 

Figures [3]|4] show that this choice of the radia- 
tive efficiency is appropriate for PBHs with masses 
Mpbh ~ 100 M0. The spectrum of the radiation emit- 
ted by the accreting material depends on the accretion 
geometry, accretion rate, and the mass of the black hole. 
For a thin disk, any given radius emits black body ra- 
diation with different emission temperature (multicolor 
disk). The disk is hot in the inner parts and colder 
in the outer parts. The maximum temperature of the 
disk depends on the black hole mass: T^ax oc M^^fj^m. 
Beyond photon energies kTmax the emission is domi- 
nated by non-thermal radiation produced by a diffuse 
hot corona around the accretion disk. The spectrum of 
the non-thermal component is a power law of the form 
i^L^ cx ly"'^ with /3 ^ 1.5. Guided by observations we 
model the spectrum as a double power-law with a break 
at energies hv = kTmax with slope (3 = 0.3 at low energy 
and /3 = 1.6 at high energy. 

Spherical accretion or ADAF: If m < 1 the accre- 
tion efficiency depends on the flow geometry. If the 
infall is quasi-spherical the radiative efficiency is min- 
imal. The dominant emission mechanism is thermal 
bremsstrahlung. Most of the radiation will originate from 
the region just outside the event horizon. For the case 
of accretion from a neutral Hi gas, including relativistic 
effects, the efficiency for con version of rest- mass energy 
into radiation is e = O.Ollrh ()Shapirolll973a lb) . Hence, 

/ — O.Ollm^ if m < 1 (spherical accretion). (30) 

The spectrum is well approximated by a power law of the 
form vLi, cx v^'^ dX v > 13.6 eV and ha s an exponential 
cut-off at - 5 X 10^ eV ()Shapirolll973al ). This case gives 
the most conservative estimate of the effect of PBHs on 
the cosmic ionization history. In addition some authors 
have found that when the effects of magnetic field and 
cosmic rays are included in the calculations, the radiative 
efficiency can be larger than the 0.01m assumed in our 
fiducial case ()Meszaroslll975b[ ). 

If the gas has non-negligible angular momentum, an 
advection dominated accretion flow (ADAF) may form 
(|Narayan &: YillTggSl ). In this case the radiative efficiency 
is a factor of ten larger: e = O.lm. Hence, I = O.lm^. 
This case is not particularly relevant for our study as 
we have seen in § 13.31 that only PBHs with masses > 
1000 Mq can form accretion disks and typically for these 
masses we find m > 1. 

3.6.1. Is the emerging radiation trapped by the accreting 

gas? 

Assuming spherical accretion and negligible pressure 
with respect to the gravitational potential, the gas falling 
onto PBHs acquires a velocity that approaches the free 
fall velocity near the black hole.^ From mass conserva- 
tion it follows that the gas density profile is pg oc r~^-^ 

^ For realistic cases in which the angular momentum of accreting 
gas is important we refer to previou s well known studies of hot 
optic ally thin accretion flows UShapiro et al.iri976t INaravan fc Yil 



in the inner parts. Hence, the gas column density 
Ng (X r^^/^ diverges for radii r ^ and so does the 
Compton scattering optical depth for outgoing radiation. 

We checked whether the accreting gas may become 
opaque to Compton scattering. In this case, the infalling 
gas would trap the emerging X-ray photons (which are 
mostly emitted near the black hole horizon) or the X-ray 
photons may be reprocessed into lower energy photons. 
The optical depth for emerging radiation emitted at a 
distance rmin — oiRsch from the black hole, 

Te^—i Pgdr = (31) 

where a > 1 is the distance in units of the black 
hole Schwarzschild radius and where we have assumed 
^max ^min- HcnCC, for any value of the accretion rate 
smaller than m ~ a^^^ ~ 1 the gas is optically thin to 
Compton scattering. Since we assume that PBHs with 
m > 1 form an thin accretion disk, we can safely neglect 
trapping of radiation emitted near the black hole in all 
the models we consider. Similar calculations show that 
the gas accreting onto PBHs is transparent (to Compton 
scattering) to external background radiation. 

4. FEEDBACK PROCESSES 

Gas accretion onto PBHs produces radiation that heats 
and ionizes the IGM, thus raising its temperature above 
the value given in equation ([5]). X-ray photons, hav- 
ing a mean free path larger than the mean distance be- 
tween the sources, tend to build up a uniform radiation 
background that increases the ionization fraction of the 
gas and heats the IGM. The UV radiation, on the other 
hand, is absorbed not very far from the emitting sources, 
so UV photons produce spheres of fully ionized hydrogen 
around each PBH. These processes introduce global and 
local feedback effects which may alter the accretion rate. 
Here we describe these effects and comment on their im- 
portance for the accretion calculation. 

4.1. Global Feedback 

We use a semianalytic c ode, described in de- 
tail in (iRicotti fc Ostrikerl l2004f) and based on 
(jChiu fc Ostrikeii 120001 ). to follow the chemical, ioniza- 
tion and thermal history of the IGM from recombination 
to the rcdshift of the formation of the first galaxies. We 
consider a gas of primordial composition and include H2 
formation/destruction processes. Although the code in- 
cludes a stellar reionization and galaxy formation formal- 
ism, in the present paper we focus on the cosmic epochs 
preceding the formation of the first galaxies. In addition 
to the emission from PBHs we include the contribution 
from the CMB radiation background that gives us a red- 
shift of recombination Zrec ~ 1000, in good agreement 
with observations. 

The thermal feedback by the X-ray background is cal- 
culated self-consistently: the temperature and ioniza- 
tion of the cosmic gas is determined by the luminosity 
of PBHs, and the PBH accretion rate is a function of 
the gas temperature and ionization. The thermal and 
chemical history is calculated iteratively until the solu- 
tion converges. The total emissivity per unit mass from 
accretion onto PBHs is proportional to I fpbh- At z < 30, 
when the age of the universe is t ~ 10^ yr, accretion onto 



10 



PBHs stops due to their peculiar velocities. Even if we 
assume that PBHs accrete at the Eddington rate before 
z = 30, the accretion time scale tsai-p = 4.4 x 10^ yr is 
much larger than the age of the universe at this redshift. 
Hence, the growth of PBHs by gas accretion is negligible. 

4.2. Local Feedback and Duty Cycle 

In addition to tracking the X-ray background, 
the semianalytic code simulates the evolution of the 
Stromgren spheres produced by UV radiation. These 
ionized bubbles have a small volume filling factor and 
thus a negligible effect on the ionization history. Never- 
theless, the gas is ionized and heated by UV radiation 
within the Hii regions, leading to local feedback effects 
that may reduce the gas accretion rate. 

Local feedback may reduce the gas accretion rate if: 
i) the radius of the Hii region exceeds the Bondi radius 
(thii > tb), and ii) if the gas temperature inside the H ii 
region is higher than the temperature outside. Although 
this is always the case at low redshifts, at z > 30 the 
Compton coupling between the free electrons inside the 
Hii regions and the CMB radiation is very effective in 
keeping the gas temperature near the CMB temperature. 

If the final size of the Hii region is larger than the 
Bondi radius {rnii/rB > 1) and the temperature inside 
the Hii region is much larger than that outside, then the 
gas accretion rate may stop or decrease. For simplicity, 
let's consider the most extreme scenario, in which when 
thii/tb ^ 1, the luminosity becomes negligible and the 
H II region recombines. With this assumption the time 
averaged luminosity can be estimated to be 



{l)t 



I 



I 



I + toff /ton l + [rHii/rBYl^ 



fdutyh (32) 



where tof / — trec,H is the H i recombination timescale 
and ton — trec,H{'fB /thiiY^^ IS thc timcscalc it takes 
for the ionization front to reach the Bondi radius r^. 
Thus, it follows toff /ton = {rnii /tbY^^ ■ Equation (l32|) 
gives a rough estimate of the duty cycle produced by 
local feedback. 

4.2.1. Temperature structure of the Hll region 

We have performed time-dependent ID radiative trans- 
fer simulations of the Hii regions around PBHs to 
estimate rnii and the temperature Thii inside the 
Stromgren radius. The basic results can be also derived 
analytically. The radius of the Stromgren sphere around 
a PBH of mass Mpf,^ at the center of a spherical gas cloud 
with density profile Ugas — nHo(r /tq)'^ , where nnoiz) is 
the mean cosmic density of the gas, is 



THII = (0.77 pc) 



5, 



0,49 



/7 



1/3 



l + z V 
1000 / 



(33) 



where = 1/3 for 7 = (i.e., assuming gas at uniform 
density) or f-^ ^ ln(ro/rmi„) ~ 10 assuming a free falling 
gas that becomes coUisionally ionized at Vmin {i-c., we 
adopt 7 = —1.5 and ro/rmin ^ 10^). We have defined 
the quantities 



-1 



(34) 



the number of ionizing photons emitted per second in 
units of 10"*^ photons s^^, and 



X = 



Lboi 



hv 



dv 



(35) 



the mean energy of the emitted photons in Rydberg units 
{hvQ = 13.6 eV). Assuming a power-law spectrum with 
index (3 = 0.5 in the frequency interval 13.6 < v < 10^ 
eV, appropriate for spherical accretion, we find X ^ 192. 
For thin disk accretion we find X = (3/ {(3 — 1) 3 for 
(3 — l.b. In order to estimate thc value of thc Stromgren 
radius we assumed an H i recombination coefficient a — 
1.28 X 10~^^ cm~'^s~^, appropriate for a gas at T ~ 2, 000 
K. This temperature is typical for the gas inside an Hii 
region at z ~ 500. 

We now estimate the temperature inside the H 11 re- 
gion assuming that the dominant gas heating is H i pho- 
toionization and the dominant cooling is Compton cool- 
ing. In a uniform density gas in ionization and ther- 
mal equilibrium the temperature inside the H II region is 
constant. At equilibrium the ionization rate per hydro- 
gen atom equals the recombination rate: tion=trec,H 
(124 yr)[(l-|-z) /lOOO]"^ and the heating rate is {hv)/t,on, 
where {tiu) ~ 0.365 eV/(/3 -I- 2) is the mean pho- 
ton energy deposited into the gas per hydrogen atom. 
Neglecting Compton heating, which is only important 
in a small volume around the black hole, a fraction 
C 0.0268/(/3 -I- 2)X of the radiation emitted by the 
black hole is deposited as heat inside the Hii region. 
If we assume that the heating rate equals the Comp- 
ton cooling rate k{THii - TcMB)/tcom,p, where tcomp ^ 
(0.76 yr)[(l + z)/1000]-4 we find 



Thii 

TcMB 



0.36 fl 



(i- 



1000 



(36) 



Thus, for spherical accretion with (3 = 0.5 (disk accretion 
with 13 = 1.5) the temperature inside the Hii regions at 
z ~ 500 is only 56% (41%) higher than the temperature 
outside, which is roughly equal to the CMB tempera- 
ture TcMB ~ 1375. We conclude that at high redshift, 
even if the Bondi radius is smaller than the Stromgren 
radius, the increase of the IGM temperature due to local 
feedback is negligible. Hence, we can neglect the local 
contribution to thermal feedback (from UV photos) with 
respect to the global thermal feedback from X-ray heat- 
ing, that is always included self-consistently (see § 14. ip 

More sophisticated radiative transfer simulations con- 
firm the results of the analytical calculations of rnii and 
Thii- In Fig. [5] we show the time evolution of the tem- 
perature profile inside an H il region around a PBH with 
50 = 1052 at z ~ 500. The ID radiative trans- 
fer code used to p roduce the tempera ture plots is de- 
scribed in detail in iRicotti et al.l (|2001f ). In addition, in 
the present simulation we have also included Compton 
heating from X-rays emitted near the PBH. The increase 
of the gas temperature at small radii evident in Fig. [5l 
T cx r~2^ is produced by Compton heating. We estimate 
that Compton heating becomes dominant over photoion- 
ization heating at radii 
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find that the ratio 
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Fig. 5. — (a) Temperature structure of the Hll region around a 
PBH at 2 = 500 emitting So = 10^2 ionizing photons per second. 
The curves show the temperature profile after 2 yr (dashed), 100 yr 
(dotted) and 4600 yr (solid) after the source turns on. The source 
spectrum is one appropriate for spherical accretion onto a black 
hole, with log-slope f3 = 0.5. 
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1000 / 



(37) 



After a time-dependent phase during which the Hii re- 
gion reaches its Stromgren radius, the temperature in- 
side the H II region decreases to a constant value Thii ^ 
2, 100 K, about 56% higher than the temperature out- 
side the Hli regions, in agreement with the analytical 
estimate. We have verified that the temperature Ty// is 
independent of the source luminosity and depends on f3 
and the redshift according to equation ((36l) . 

4.2.2. Stromgren vs Bondi radii 

Let's now estimate the ratio rjiji/rg, which, along 
with the HII region gas temperature, determines whether 
local feedback will reduce the gas accretion rate. ^From 
equation ^ and equation psp with /-y = 10 we find 



TB 



'3x 10" ( - 



l + z 



M, 



pbh 



1 Mr, 



1000 



5.7km 



(38) 



for "naked" PBHs. If we include the growth of an ex- 
tended dark halo, rnii /tb decreases by a factor 0.3[(1 -I- 
z)/1000] at z < 1000. We shall consider two cases 
for the accretion luminosity: Case A, appropriate for 
small accretion rates m < 1 and Case B, appropriate 
form - 1-10 (sec 

Case A: I oi fn? . 
The case / cx rn? is the most interesting because is 
the most likely to occur for PBHs and because the ra- 
tio rHii/^B is independent of the gas effective velocity 



and PBH mass. Indeed, run oc /^/^ 



cx Wg^j and 



Assuming e — 0.01 (spherical accretion) we 



TB 



4- 

V200 



(39) 



is a constant of order of unity for a "naked" PBH with 
mass Mpfc^ < M„ - 10'' M© and m < 1. For PBHs 
more massive than M^r, the constant is smaller by a 
factor {Mpbh/Mcr)'^^'^ , thus local radiative feedback can 
be safely neglected because run < tb- Including the 
growth of the dark halo around PBHs we find: 



< 1, 



(40) 



valid only at redshifts z ^ 100. At lower redshifts, 
THii/rB increases with respect to the estimate in equa- 
tion (|40|) but its value remains on order unity or smaller 
for m < 1. 

Case B: I cx m. 
This case applies to the case of relatively massive PBHs 
with typical masses > 10^ that have 1 < m < 10 or 
0.1 < ^ < 1. It is sufficient to consider the case of a PBH 
clothed in its dark halo. We find 



THII 
TB 



V lO^Mp 



-2/3 



-1/3 



(41) 



Thus, rnii < tb for Mpf,/i > 10^^ Mq. Also in this 
case we can neglect local feedback with the exception 
of the mass range 10^ — 10^ M© where local thermal 
feedback may reduce accretion by a factor of a few (from 
equation (j32p we obtain a reduction fduty ^ 30% for 
Mpbh ~ 10^ Mq). 

In summary, the local feedback due to the formation of 
an H II region around the PBH can be neglected in most 
cases. At high redshift the temperature inside the Hii 
regions is only a few tens of a percent higher than the 
temperature outside (see equation [35]), thus even if the 
Bondi radius is smaller than the Stromgren radius, the 
reduction of the accretion rate due to thermal feedback is 
negligible. At lower redshifts [z < 100) the temperature 
inside H ii regions is larger than the temperature out- 
side but we found that the Stromgren radius is typically 
smaller than the Bondi radius. 

The only exceptions are massive PBHs with IO^Mq < 
Mpbh < 10^ M0 for which we estimate fduty ^ 30%. For 
this case, since thii > tb, we need to also assume that 
the gas is fully ionized at the Bondi radius. Hence, in ad- 
dition to thermal feedback, the ionized gas inside the H II 
region increases the gas viscosity due to Compton drag. 
Compton drag is negligible at z < 100 (see Figs. [3]|4]) 
but reduces the accretion rate at z > 100 when, instead, 
thermal feedback is negligible. 

5. SIMULATING THE COSMIC IONIZATION 

We simulate the ionization, chemical and thermal his- 
tory of the universe after recombination using a modi- 
fication of the semianaly tic code discussed in detail in 
iRicotti fc OstrikeH (pOOl . The original code is tailored 
to simulate UV and X-ray reionization from high redshift 
galaxies at z < 30. Here we do not include any ionizing 
source other than PBHs. In order to simulate standard 
recombination properly we model the redshifted CMB 
black body radiation including the effect of Thompson 
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Fig. 6. — Le/i. Simulation of early ionization by PBHs with mass Mpi,^ = 100 Mq and abundance = 10 The panels show: (a) 
the ionization fraction Xe (the dashed line shows the recombination history without PBHs); (b) the temperature (the dashed line shows 
the thermal history without PBHs); (c) the He ionization fraction; (d) the H2 abundance; (e) the accretion rate, m (solid curve) and 
luminosity / (dashed curve); (f) the optical depth to Thompson scattering t^. Right. Same as in the left panel but for Mpi,/^ = 1000 Mq 
and fpijh = 10"'' (solid curves) and /pjh = 10~^ (dotted curves). 

ton cooling and molecular hydrogen cooling. The vol- 
ume filling factor of the coUisionally ionized gas is fv ~ 
npbhiTB , where ripb/i is the PBH physical number 
density. Assuming that a fraction of the IGM fy is 
fully ionized and the rest of the cosmic gas is neutral 
we have a volume-weighted ionization fraction due to 
shocks: {xe)v ~ lO-^^fpbh{Mpbh/l Mq)\ This effect 
is negligible for any PBH mass assuming values of fpbh 
allowed by observational constraints. 

Examples of simulations with Mpbh = 100 Mq and 
1000 Mq are shown in Fig. [6l Here, we describe sim- 
ulations in models that are consistent with FIRAS and 
WMAP3 data. In these models the contribution of PBHs 
to the total optical depth to Thompson scattering typ- 
ically does not exceed Are ~ 0.03 — 0.04 at 68% confi- 
dence level and Are ~ 0.07 — 0.09 at 95%confidence level. 
The contribution is uncorrelated with the value of Te pro- 
duced by ionization sources in high redshift galaxies (at 
z < 30) which is re = 0.09 ± 0.03 (see § E]). 

The mean electron fraction increases approximately as 
Xe oc {z + from Xe ^ 10^^ at 2: '--^ 900 to values close 
to Xe ^ 10^^ — 10^^ at z ^ 10. Hence, the viscous effect 
of Compton drag can be neglected after recombination. 
The main contribution to the partial ionization of the 
cosmic gas after recombination is from X-ray emission at 
redshifts z > 100. At smaller redshifts, due to the rapid 
decline of the accretion rate onto PBHs, the ionization 
fraction keeps increas ing slowly, due to r edshifted X-ray 
background photons (jRicotti fc Ostrikei1 [2004). The gas 
temperature is approximately equal to the CMB temper- 
ature at z > 200 and ranges between 100 K to 1000 K 
afterwards. The He 11 ionization fraction becomes 1 — 10% 
and the molecular abundance XH2 ~ 10"'* — 10~^ after 
redshift z ^ 100. The molecular abundance is 10-100 
times larger than the standard value, XH2 ~ 10~^, ob- 



opacity. We test the recombination history calculation 
in the absence of PBHs ag ainst a widely used code REC- 
FAST (jSeager et al.lll999i ). We can reproduce the recom- 
bination history at redshifts z < 900 (when < 0.1) 
with accuracy of a few percent. Using the equations de- 
rived in this work we model the UV and X-ray emission 
from accreting PBHs including feedback effects, the sec- 
ondary ionization due to fast photo-electrons, and Comp- 
ton heating/cooling, and we solve the chemical network 
for a gas of primordial composition {e.g., molecular hy- 
drogen formation/dissociation). We run a grid of models 
with a range of PBH masses spanning several orders of 
magnitude, between 10"^ Mq to 10^ Mq. 

For the sake of simplicity, throughout the rest of the 
paper we consider the effects of a population of equal 
mass PBHs which account for a fraction fpbh of the dark 
matter. The results can be used to estimate the effect 
of an arbitrary PBH mass function by integrating the 
optical depth dre (M) in each mass bin with weight func- 
tion {d/pbh/dMy/"^. See below for an explanation of the 
choice of the weighting function. Special care should be 
taken in excluding from the integration PBHs which are 
sufficiently small to be part of the dark halos of more 
massive PBHs. In addition we need to check that that 
the proper velocities of small PBHs, due to their interac- 
tion with more massive ones, is negligible (i.e., subsonic). 
Fitting formulas for as a function of Mpbh and fpbh are 
provided at the end of this section. 

Independently of the radiation emitted by the gas 
falling directly onto the PBHs, gas inflow onto the PBH 
can produce shocks and coUisional ionization. For the 
idealized case of an isothermal gas the sonic point is 
located at O.Srs, where is the Bondi radius. At 
this point a shock develops that ionizes the gas, even if 
the temperature remains below 10,000 K due to Comp- 
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tained neglecting PBHs. Depending on the mass of PBHs 
we encounter different regimes for the accretion. 

1) 1 Mq < Mpbh < 30 M0 : The dark halo has negligible 
effect on the accretion rate. The accretion is spherical 
and the luminosity / (x m? . Global feedback effects are 
negligible. 

2) 30 M© < Mpbh < 300 M© : The dark halo increases 
the accretion rate onto PBHs, mostly at redshifts z < 
100. The accretion is spherical and the luminosity / cx 
rn^. Global feedback effects are important in reducing 
the accretion rate aX z < 100. 

3) 300 M© < Mpbh < 3000 M©: The dark halo increases 
the accretion rate onto PBHs at all redshifts. The ac- 
cretion rate is large (m ^1), thus a thin disk may form 
and the luminosity is / cx fdutyfn- Global feedback effects 
moderately reduce the accretion rate. 

4) Mpbh > 3000 M© : The dark halo increases the 
accretion rate onto PBHs at all redshifts. For masses 
> 5 X 10^ M©, the Hubble expansion becomes important 
and the Bondi solutions transition to self-similar infall 
solutions (see Paper II). PBHs accrete at the Edding- 
ton rate with a duty cycle fduty Global feedback effects 
do not reduce the accretion rate significantly for values 
Are < 0.05, compatible with WMAP3 (see §0- 

The qualitative dependence of Xi,{z) and Te(z), on the 
mass and abundance of PBHs can be understood assum- 
ing ionization equilibrium. For Xe{z) <gC 1 and neglecting 
the temperature dependence of the recombination coef- 
ficient we have Xe{z)'^ cx I {z, Mpbh) f pbh- For PBHs with 
masses smaller than Mpbh ~ 100 M© we have m < 1 and 

/ cx Mpjj,^. Hence, Te cx Xe cx Mpbhfpl^- If the mass of the 
PBH is ^ 100 M© the accretion rate is typically rh > 1, 
I cx Mpbh and Tg cx a;e oc (Mpbhfdutyfpbh)^^'^ ■ Finally, if 
Mpbh ^ 1000 M© we have rh > 10, PBHs accrete nearly 
at the Eddington limit and Tg cx Xg cx {fdutyfpbhY^'^ is 
independent of the PBH mass. To first order, the value 
of Are produced by PBHs of mass Mpbh and abundance 
fpbh can be parameterized as follows: 



0.05 



Are 



(^)4. ifMp,.<100M©, 
0-1 (^) ifl02<i^<103 
I IQ^ {fduty fpbh)'^ if Mpbh > 10^ M, 



0- 

(42) 

The fit is derived from a grid of simulations with Are < 
0.2 and is not accurate for large values of Are due to 
feedback and saturation effects. 

6. EFFECTS OF PBHS ON THE CMB SPECTRUM AND 
ANISOTROPIES 

6.1. Temperature and Polarization Anisotropies 

In § [S]we have shown that a signature of the existence 
of PBHs is the modification of the cosmic recombination 
history. It is possible to construct models with re ~ 1 due 
to the large residual fractional ionization of the cosmic 
gas after recombination produced by X-ray ionization. 
Of course, such models are ruled out by CMB observa- 
tions. The WMAP3 constraint on is Te ~ 0.09 ± 0.03. 
This limit on assumes that the intergalactic medium 
becomes partially or fully ionized by stars and/or black 
holes starting at redshift z ^ 30. As illustrated in Fig. [71 
the partial ionization from PBHs and the ionization from 



high redshift galaxies forming at z < 30 produce very dif- 
ferent signatures on the CMB polarization anisotropies. 
PBHs affect small angular scales with I ^ 10, while radi- 
ation emitted by high redshift galaxies affect larger scales 
with I < 10. Hence, PBHs do not contribute significantly 
to the value re = 0.09 quoted by the WMAP3 team, 
which must be produced by ionization sources other than 
PBHs. 

In order to constrain models which allow for the ex- 
istence of PBHs we include an additional cosmologi- 
cal parameter describing the deviation from the stan- 
dard recombinatio n history calculated using RECFAST 
(jSeager et al.ll999l) . We have modified th e publicly avail- 
able codes CAMCMB and COSMOMC (iLewis fc Bridld 
|2002f) to include an additional degree of freedom. The 
redshift dependence of the ionization history after re- 
combination depends weakly on the mass of PBHs for 
values of Are ~ 0.1. We modify the recombination his- 
tory Xe.reciz) givcn by RECFAST as follows: 



Xe{z) 



XeO 



l + z \ 
1000 ) 



,0.1 



(43) 



The free parameter Xeo is the ionization fraction at 
z — 1000. The ionization fraction increases proportion- 
ally to the scale parameter after recombination and is 
constant at later times. In addition, we model the in- 
stantaneous reionization at 2 = z^ei produced by galactic 
sources. We will show that the effect on the CMB polar- 
ization anisotropies of PBHs and other ionizing sources 
are uncorrelated. 

Using the Markhov Chain Monte-Carlo code COS- 
MOMC with WMAP3 data we find Icr and 2a marginal- 
ized confidence limits for the new free parameter: Xeo < 
2 X 10-"^ (68% CF) and Xeo < 4.2 x 10"^ (95% CF). 
These upper limits on the fractional ionization corre- 
spond to upper limits on Are = Te — Te rei of 0.05 (68% 
CF) and 0.1 (95% CF). As illustrated in Fig.li the new 
parameter does not correlate with the redshift of reion- 
ization Zrei but because it increases the total it corre- 
lates with (Tg and ris , increasing their value s with respect 
to th e ones quoted by the WMAP3 team (iSpergel et al.1 
I2006f) . Invoking a non-standard recombination history 
may ease the tension between the low value of erg ~ 0.74 
from the WMAP3 analysis (assuming standard recom- 
bination history) and clusters d ata that in stead seem to 
favo r larger values of erg ^ 0.9 (lEvrard et al. 2007). but 
see (|Bode et al.ll2007( ). 

Using the results of the simulations presented in § O 
that are approximately summarized by equation (|42l) , we 
are able to constrain the abundance of PBHs with masses 
1 M© < Mpbh < 10® M©. The results are shown by the 
thick solid line in Fig. [9fleft). 

6.2. Spectral Distortions 

The injection of some form of energy density, A[/, 
in the expanding universe at redshifts z < 10^ can- 
not be fully thermalized, and so it produces observ- 
able deviations of the CMB spectrum from a perfect 
black-body (e.g.. iBurigana et al.lll99lh . There are two 
types of spectral distortions: ^ and y-distortions. En- 
ergy injection in the redshift range 10^ < z < 10^ pro- 
duces a Bose-Einstein spectrum with chemical potential 
AU/U — 0.71/i, where U is the unperturbed energy den- 
sity of the CMB. In the redshift interval 10^ < z < 1000 
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Fig. 7. — Power spectrum of temperature and polarization 
anisotropics for the best fit WAMP3 model witli = 0.09 and 
redshift of reionization z^ei = 11 (crosses) and a model with the 
same Te but z^ei = 7 and modified recombination history with 
constant residual ionization fraction x^o = 2.5 X 10""^. The panels 
show from top to bottom: TT, EE and TE power spectra. 

the Coniptonization distortion dominates and has am- 
phtude AC//C/ = Ay. Using the spectrom eter FIRAS on 
the COBE satelhte, iFixsen et all ()1996D found the fol- 
lowing upper limits for the deviation of the CMB from 
a Planck spectrum: ^ < 9 x 10~^ and y < 1.5 x 10~^ 
at at 95% confidence. It follows that MJ/U < 6 x 10"^ 
at 95% confidence from 10^ < z < 1000. In this section 
we calculate the value of the y-parameter produced by 
a population of PBHs of mass Mp^h and compare it to 
observed upper limits from FIRAS to set upper limits 
on their abundance. In order to calculate the total value 
of the y-parameter we need to distinguish between three 
epochs for the energy injection: 

1) Before the redshift of last scattering at Zrec ~ 1000 the 
universe is optically thick to Compton scattering hence 
all the energy emitted by accreting PBHs is absorbed by 
the cosmic gas. During this epoch the y-parameter is 

^ ^ dz dAUjz) 

" 4C/(zeg) J,^^ aH{z) dt ^ ' 

where dAU{z)/dt is the total energy per unit comoving 
volume per unit time emitted by accretion onto PBHs 
and U{z) is the unperturbed comoving energy density of 
the CMB at redshift z. We will show that this contribu- 
tion to the total y-parameter is dominant. We integrate 
dAU(z)/dt starting at the redshift of matter-radiation 
equality because before z^q the effective cosmic sound 
speed approaches the value Cg ~ and the gas ac- 

cretion rate onto PBHs is negligible. 

2) Between last scattering and decoupling the estimate 
of the y-parameter is complicated by the fact that only 
a fraction of the energy emitted by PBHs is absorbed by 
the gas and exchanged with the CMB radiation. As a 
zero-th order approximation we can use equation (|44p . 
integrating between z^ec and Zrec, to estimate an upper 



limit for the y-distortion. The integration shows that 
for PBH masses Mpth < 10"^ Mq, y2 <C yi even when 
we assume that all the energy is absorbed by the gas. In 
i^ l4.2.1l we estimated that only a fraction C ^ 1.4x 10^** of 
the radiation emitted by PBHs is deposited into heat. A 
more realistic estimate of AU is obtained by multiplying 
our upper limit for 1/2 by C. Hence, we do not need to 
worry further about getting a more precise estimate of 
y2 for PBH of any mass as y2 <S^yi- 
3) After decoupling Compton heating/cooling becomes 
negligible and we can calculate the y-parameter using 
the relationship: 

ys = n / [Te - Tcmb)dTe. (45) 

Let's estimate the upper limit for y^. Assuming 
Tf, = constant ^ Tcmb we have ys = 1.726 x 
10-^(re/10'' K)Are. In all the models consistent with 
WMAP3 data the gas temperature is Te < 1000 K and 
Are between redshift z ~ 10 and z ~ 100 is < 0.05. It 
follows that the upper limit on ya is 8.6 x 10~^ which 
is negligible when compared to the FIRAS upper limit. 
Finally, the value of the y-parameter is the sum of the y 
during each epoch. Thus, in our case, we have y ~ yi- 

Let's now estimate yi. The total energy density emit- 
ted per unit time per unit comoving volume is AU /dt — 
ILEdJipbh- Thus we find 

AMp,„anT,fil + Zeq)^P'''J^^^ aH{z) ' 

(46) 

Using equation ([50]) for the dimensionless accretion lu- 
minosity I (see also Fig. S]), we obtain the value of the y- 
parameter as a function of Mp^h and fpbh ■ Imposing y < 
1.5 X 10^^ we obtain upper limits for fpbh{Mpbh) at 95% 
confidence. The results are summarized in Fig. [5Kleft). 

In summary, before the redshift of recombination gas 
accretion onto PBHs with mass < 100 Mq is not greatly 
reduced by Compton drag. Although the accretion lumi- 
nosity during this epoch does not contribute to increase 
Te, the energy injection produce spectral distortions of 
the CMB, increasing the value of the y-parameter. We 
find that the existence of PBHs with masses < 100 Mq 
is best constrained by upper limits on the y-parameter 
from FIRAS. 

7. DISCUSSION AND SUMMARY 

During the radiation era, mildly non-linear perturba- 
tions with 6p/ p ^ 0.5 — 1 can collapse directly into pri- 
mordial black holes (PBHs). Such black holes may have 
masses ranging from the Planck mass to a million solar 
masses, depending on the redshift of their formation and 
the details of the formation mechanism. The abundance 
of evaporating PBHs with masses < 10""^^ g is constrained 
by observations to be a fraction (3 ^ 10~^^ of the mean 
energy density of the universe at the time of their forma- 
tion, but the existence of PBHs with masses larger than 
10^^ g is poorly constrained. It is not ruled out that the 
bulk of the dark matter may be composed of PBHs with 
masses in the range between 10^*^ g and 10^^ g or Planck 
mass relics with M ~ 10~^ g. In this work, the last of 
a series of three papers, we study the effects of a (yet 
undetected) population of non-evaporating PBHs on the 
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Fig. 8. — Likelihood isocontours of cosmological parameters as a function of x^q which describes the modified recomb ination history 
produced by PBHs (see the text). We have used WMAP3 data and a modified version of COSMOMC I ILewis fc Bridlell200l ) . The equation 
used for the non-standard recombination history is given in equation 1 143 I I and the free parameter x^o describes the electron fraction at 
z = 1000. The four panels from the top-left corner show (clockwise): i) marginalized likelihood of XeO (in units of 10~*) with 68% and 
95% confidence limits; ii) optical depth Te vs Xeo; iii) amplitude of the power spectrum vs Xeo; iv) scalar tilt of the initial power spectrum 
depth Us vs XeO- 




thermal and ionization history of the universe and their 
signatures on the CMB anisotropies and spectrum. In 
Paper I we focused on studying the formation and growth 
of the dark matter halo which envelopes PBHs that do 
not co nstitute the b ulk of the dark matter. In the second 
paper (jRicottil [20071 ). we study in detail the Bondi-type 
accretion solutions onto PBHs including the effects of 
Compton drag, Hubble expansion and the growth of the 
dark matter halo. Finally, this work focuses on model- 
ing the accretion luminosity of PBHs including feedback 
effects and observational signatures. 

We find that if a fraction fp^h of the dark matter 
is in PBHs with mass > 0.1 Mq, the energy released 
due to gas accretion may produce spectral distortions 
of the CMB radiation and keep the universe partially 
ionized after recombination. The limits on the mass 
and abundances of PBHs set from observations of the 
X-ray background are much less restrictive than those 
from the CMB. The modified recombination history pro- 
duces observable signatures on the spectrum of polariza- 
tion anisotropies of the CMB at angular scales I ^ 10. 
Hence, the effect of PBHs cannot be confused with the 
effect of ionization by high redshift galaxies which affect 
polarization anisotropies on larger angular scales. 

We are able to improve the constraints on fpbh for 
PBHs with masses > 0.1 M© by several orders of mag- 
nitude using WMAP3 and FIRAS data. The results are 
summarized in Fig. [5] (left). The upper limits on the 
abundance of PBHs with masses O.IM© < Mpth < 10^ 
Mq at the epoch of their formation, /3, are shown in Fig. [5] 



(right). We use equation ([T]) to derive /3 as a function of 
Mpbh, fpbh and the ratio fnor = Mpbh/Mn between the 
PBH mass and the mass Mh of the horizon at the epoch 
of PBH formation. 

Fitting WMAP3 data with cosmological models that 
do not allow for non-standard recombination histories as 
produced by PBHs or other early energy sources may 
lead to an underestimate of the best-fit values of the am- 
plitude of linear density fluctuations, cg, and the scalar 
spectral index, n^. This happens because the contribu- 
tion of PBHs to the optical depth to Thompson scat- 
tering, which is uncorrelated with the contribution from 
galactic ionization sources, can be Are ~ 0.05 (Are ^ 0.1 
at 95% CF). Since rig and are correlated with Te, their 
best fit values may increase to ~ 1 and erg ~ 0.9 (at 
95% CF). This is a general result that may reduce recent 
tensions bet ween WMAP3 data and clusters data on the 
value of as ()Evrard et al.ll2007D . 

A population of intermediate mass black holes 
(IMBHs) with masses of about 100 - 1000 Mq is still 
allowed and may be widespread if a fraction of the ul- 
traluminous X-ray s ources (ULX) observed in nearby 
galax ies host IMBHs (jMiller et al.ll2003HDewangan et all 
I2006D . The origin of IMBHs is unknown, but if they 
are produced by Pop III stars their number may fall 
short in explaining t h e observed ULXs popu lation {e.g., 
iKuranov et all l2007t iPelupessv eFall |2007[ ). We find 
that, if all or a fraction of observed ULXs are PBHs with 
masses Mpbh ~ 100 - 1000 Mq with fpbh ~ 10-^ they 
would increase the best fit value of the optical depth to 
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Fig. 9. — (left) Upper limits on the present abundance of PBHs. The thick lines are the results obtained in the present work. The solid 
lines show the upper limits using WMAP3 data (CMB anisotropics) for two values of the black hole duty cycle fduty = 1 and 0.1. The 
dashed lines show the limits using FIRAS data (CMB spectral distortions) at 95% and 68% confidence. The other lines refer to previous 
upper limits from microlensing (EROS and MACHO experiments) and dynamical constraints (see introduction). (Right) Upper limits on 
the abundance of PBHs at the epoch of their formation /3 as a function of their mass. We assume that the mass of PBHs is a fraction fnor 
of the mass of the horizon at the epoch of their formation. The thick curves show the upper limits obtained in the present work and the 
thin dotted curve are limits from the EROS collaboration (microlensing experiment). 



Thompson scattering to Tg ~ 0.2. Since the scalar spec- 
tral index and the amphtude of density fluctuations 
As and erg are correlated to Tg, their best fits also increase 
to ~ 1 and as ~ 0.9. PBHs in this mass range may be 
produced in two-stage inflationary models designed to fit 
the low WMAP quadrupole (Kawasaki et al. 2006). Wc 
emphasize again that this effect is more general than the 
specific case of PBHs discussed in this paper. Any mech- 
anism or energy source that modifies the standard recom- 
bination history may affect the estimate of cosmological 
parameters in a way similar to that discussed here. 

Our results are in contradiction with the suggestion 
that MAC HOs are PBHs with mass 0.1 - 1 M© and 
fpbh ^ 0.2 (jAlcock et al.ll2000D . Such a PBH population 
would produce spectral distortions incompatible with FI- 
RAS data. 

The luminous QSOs found by SLOAN at z '--^ 6 are 
thought to be powered by 10^ - 10^ Mq SMBHs. It 
is difficult to produce such massive black holes starting 
from small seeds by gas accretion because the age of the 
universe at 2; = 6 is a few tens the Salpeter accretion 
timescale. A few massive PBHs or numerous less mas- 
sive PBHs may help explain the origin of SMBHs at high 
redshift and in present day galaxies by producing rela- 
tively massive "seeds" . Are the upper limits on the num- 
ber of PBHs derived in this work compatible with this 
scenario? The fraction of mass in SMBHs today is ap- 
proximately ^Ismbh/^dm 2.13 X 10^^ (Gcbhard t et al.l 
120001 iRicotti fc Ostrikeii [200l . For PBHs with mass 

> 1000 Mq we found fpth = %bhl^dm ^ /fduty ■ 

Hence, assuming that only a fraction Fagn < 1 of PBHs 
is incorporated into SMBHs and grows by gas accre- 
tion by a factor Xacc > 1 we have: fpbhXaccFagn ~ 
2 X 10"^ or XaccFag7i ^ '^Ofduty The most massive 
PBHs have Fagn 1 because they spiral in to the 
centers of galaxies by dynamical friction on a shorter 
timescale (tfrtc/tH{z) ^ Q.Q2Mhaio{z) / Mpbh, where Ih 
is the Hubble time) and because they may accrete gas 
more efficiently. Hence, for fduty ~ 3% and Fagn — 1 
we find Xacc ^ 1 indicating that even scenarios with 
negligible mass accretion onto PBHs {i.e., only growth 
through mergers) are consistent with the observed mass 



in SMBHs today. 

Less massive PBHs have lower probability for growing 
to masses typical of SMBHs because the Bondi accretion 
rate is oc M^. However, the upper limit on the abun- 
dance of PBHs increases steeply with decreasing mass 
for Mpbh < 1000 Mq. Thus, although a smaller fraction 
of the seed PBHs can grow substantially, the number of 
seeds available can be much larger. PBHs with masses 
smaller than 100 Mq, assuming Bondi type accretion 
from the ISM of a typical high-z galaxy, are unlikely to 
accrete rapidly enough to grow to SMBH masses in less 
than 1 Gyr, even if they constitute a few per cent of the 
dark matter dKuranov et al.l[2007l : iPelupessv et al.ll2007l : 
IRicotti fc Kockertj|2007[ )" 

The increased fractional ionization of the cosmic gas 
produced by non-standard recombination also increases 
the primordial molecular hydrogen abundance to XH2 ~ 
10"'' - 10~5 after redshift z ~ 100. This value is be- 
tween ten and one hundred times larger than the stan- 
dard value, XH2 ~ 10~^, obtained neglecting PBHs. The 
increase of the cosmic Jeans mass due to X-ray heating 
is negligible for models consistent with the CMB data. 
Therefore, the formation rate of the first galaxies and 
stars may be enhanced if a population of PBHs exists. 
Several aspects of first-star and galaxy- formation physics 
would be affected by the enhanced molecular fraction: 

(i) the mass of the first stars may be re duced due to 
formation of HD molecules ( Nagakura fc Om ukai 2005) ; 

(ii) the intergalactic medium would be optically thick 
to H2 photo-dissociating radiation in the Lyman- Werner 
bands, allowing molecular hydrogen to survive in the low 
density IGM even at relatively low redshifts z ~ 10 — 15; 

(iii) the epoch of domination of the first stars and galax- 
ies would probably start earlier and perhaps last longer. 
The number of first galaxies that remain completely dark 
would be reduced. It is not obvious that the star forma- 
tion efficiency and other internal properties of the first 
galaxies would be affected because feedback effects such 
as photo-evaporation from inter nal sources and SN ex- 
plosions are probably dominant (jRicotti et al.ll2002al |bh. 
We leave quantitative calculations on the impact of PBHs 
on the formation of the first galaxies to a future work. 
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